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CN \ ABSTRACT 

We investigate the properties of the pairwise velocity dispersion as a function of galaxy 
' luminosity in the context of a halo model. We derive the distribution of velocities of 

1^ \ pairs at a given separation taking into account both one-halo and two-halo contribu- 

tions. We show that pairwise velocity distribution in real space is a complicated mix- 
ture of host-satellite, satellite-satellite and two-halo pairs. The peak value is reached at 
around l/i~^Mpc and does not reflect the velocity dispersion of a typical halo hosting 
f-*) ■ these galaxies, but is instead dominated by the satellite-satellite pairs in high mass 

\ clusters. This is true even for cross-correlations between bins separated in luminos- 

. ity. As a consequence the velocity dispersion at a given separation can decrease with 

luminosity, even if the underlying typical halo host mass is increasing, in agreement 
with recent observations. We compare our findings to numerical simulations and find a 
good agreement. Numerical simulations also suggest a luminosity dependent velocity 
O ' bias, which depends on the subhalo mass. We develop models of the auto- and cross- 

correlation function of luminosity subsamples of galaxies in the observable r-pi-oj — tt 
^2 ' space and calculate the inferred velocity dispersion as a function of wave vector if 



o 



dispersion model is fit to the redshift space power spectrum. We find that so derived 
pairwise velocity dispersion also exhibits a bump at A: 1/i/Mpc. 



' Key words: cosmology: large scale structure of Universe 



1 INTRODUCTION 

Large redshift surveys such as SDSS (York et al. 2000) and 
2dF (Percival et al. 2001) have provided us with a unique 
tool to probe the 3-dimensional structure of our Universe. 
The radial distances in those surveys are, however, mod- 
ulated by the components of peculiar velocities along the 
line of sight. On large scales these velocities can be mod- 
elled using linear theory, as developed by Kaiser (1987). On 
small scales the velocities lead to radial stretching of groups 
and clusters, the so-called "fingers-of-God" (FOG). In prin- 
ciple this anisotropy in the observed structure can give us 
information on velocity structure of the Universe in addi- 
tion to the spatial distribution of objects. Peculiar velocities 
of individual objects are sensitive to the amount and distri- 
bution of dark matter, but are difficult to measure in the 
absence of reliable distance estimators. The projected pair- 
wise velocity differences, on the other hand, are much easier 
to measure from those surveys. Understanding the pairwise 
velocity distribution on all scales could therefore provide im- 
portant dynamical information on the relation between the 
galaxies and surrounding dark matter. 



The simplest approach to use dynamical information 
from pairwise velocities is to select isolated halo systems by 
some scheme (McKay et al. 2002; Prada et al. 2003; Brain- 
erd & Specian 2003; van den Bosch et al. 2004). This process 
usually results in a selection of several system composed of 
a host galaxy and their satellites. Models can then be tested 
when many such isolated systems are compared "on aver- 
age" to the theoretical predictions. The main advantage of 
this approach is that the selection scheme already labels the 
host and satellite galaxies, making the modelling somewhat 
easier. However, this approach also has several disadvan- 
tages: the number of suitable galaxies is very small (of the 
order of 1% of the total number of galaxies), the treatment 
of interlopers (galaxies that appear to be part of the system 
due to projection effects) can be diflicult and model depen- 
dent, some of host galaxies may themselves be satellites of a 
larger system and, lastly, selecting isolated galaxy-satellite 
systems may not be representative of the dark matter-galaxy 
relation in general. 

The alternative approach is to do as little preselection 
as possible and instead model the effects statistically within 
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a suitable model. This has the advantage of having large 
statistics and being more representative, given that by defi- 
nition all galaxies are being used. This is the approach taken 
in this paper. We focus on the small scale correlations in red- 
shift space and use the halo model (see e.g., Cooray & Shctli 
2002) to provide an interpretation of the observations in 
the context of a physical model. Our approach thus focuses 
on statistical properties of the entire distribution of galax- 
ies. Since the standard host-satellite separation is based on 
luminosity, this suggests that cross-correlation analysis be- 
tween faint and bright galaxies might be especially effective 
in selecting central-satellite pairs and providing information 
about the halo structure. Wc thus follow Guzik & Scljak 
(2002) in working with narrow luminosity bins and perform 
both auto and cross-correlation analysis between them. 

We assume the galaxies are either central galaxies in a 
halo or satellites, both of which need to have a specified con- 
ditional halo mass probability distribution. In this model the 
contributions to pairwise velocities come from host-satellite 
(host-sat from now on) and satellite-satellite (sat-sat from 
now on) pairs within the same halo and from two halo pairs 
(both host and satellite) . A related approach with luminosity 
thresholds has been developed by Zheng ot al. (2004), while 
van den Bosch et al. (2004) have focused on conditional lu- 
minosity function as a function of halo mass. We first build 
a model to describe the velocity structure of galaxies in the 
dark matter halos (Section 2). We perform this in terms of 
velocity distribution function, the probability that a pair of 
galaxies belonging to a given luminosity bin can be found 
at a given separation and has a given velocity component. 
This procedure is in spirit similar to that of Sheth (1996), 
although with a different emphasis and a more up to date 
version of the halo model. Wc show that less luminous galax- 
ies can appear to move faster than their typical parent halo 
population even in the simplest halo model, in agreement 
with some observational results (Jing & Borncr 2004). 

In Section 3 we compare our model with numerical sim- 
ulation. We calculate the correlation function in the rproj — tt 
plane and use it to predict the dependence of ai2 parame- 
ter one would obtain by fitting it with a simple dispersion 
model (Section 4). The last section summarises our find- 
ings and concludes the paper. Throughout this paper, we 
assume a concordant fiat cosmology with Qm = 0.3, Hubble 
constant of 70km s~^Mpc~^ and as = 0.9. 



luminosity bins of unit magnitude size. We always refer to a 
luminosity bin by its faint end magnitude. Additionally, we 
will refer to a sample of galaxies selected so that each pair 
contains one galaxy from luminosity bin Li and one galaxy 
from the luminosity bin Lj as (Li, Lj). We will use symbols 
nc{M\Li) and ni<iciM\Li) to denote the expectation value 
of the number of galaxies per halo which belong to the lumi- 
nosity bin Li and occupy central and non-central positions, 
respectively. Since every halo is assumed to have a central 
galaxy, the following must be true: 



J2nc{M\Li) = l, 



(2) 



although some of these galaxies will fall outside the lumi- 
nosity range we consider here. 

Wc assume that the central galaxy is at the centre and 
at rest with respect to the parent halo. The satellite galaxies 
are distributed in the parent halo around the central galaxy. 
We assume this radial distribution follows the NFW profile 
(Navarro et al. 1996). We denote the probability of finding a 
galaxy at a distance between r and r -|- dr in a halo of mass 
M as $i: 



$i(r|M) oc — ( H- — 



(3) 



where Ts = rvir/c, where c is concentration and rvir(Af) is 
virial radius, which is defined as a radius of a sphere of mass 
M which has the average density 200 times critical density 
of the Universe. In simulations the subhalo distribution is 
less concentrated than the dark matter. Whether the same 
is true in the real universe is still a matter of an ongoing 
debate (Gliigna et al. 2000; Berlind & Weinberg 2001; Dic- 
mand et al. 2004; CoUister & Lahav 2005; Yang et al. 2005; 
Nagai & Kravtsov 2005). In this work we assume the concen- 
tration parameter for the satellite galaxy distribution Cg = 3 
and a cut-off at the rcut = 1.5rvir, which is in a reasonable 
agreement with simulations. Simulations show a weak de- 
pendence of the best-fit concentration on luminosity, but we 
neglect this effect here. 

The probability of finding a galaxy belonging to the 
luminosity bin Li in a halo with mass between M and M + 
dM can be split into contribution from central and satellite 
galaxies: 



2 PAIRWISE VELOCITIES IN HALO MODEL 

We follow the standard halo model, which postulates that all 
galaxies live in the dark matter halos. The mass function r] = 
An /AM of dark matter halos is most reliably estimated from 
numerical simulations and here we use the mass function fit 
from simulations (Jenkins et al. 2001) 



An 



3.8-, dlno- ^ 
^ d log M ' 



(1) 

where a is the variance in the linear density field at a 
given redshift, after smoothing with a spherical top-hat fil- 
ter which encloses mass M in the moan. Each halo has a 
central gala^xy and a Poisson distributed number of satel- 
lite galaxies. We consider galaxies with r-band luminosity 
between L = —19 and L = —22 and split them into three 



P{M\Li)AM = (1 - a)PciM\Li) + aPNc{M\Li), (4) 

where a is the fraction of satellites among galaxies of lumi- 
nosity Li. 

The probability distribution for the central galaxy is of- 
ten assumed to be a 5 function. Here we assume a slightly 
more realistic model and introduce a Gaussian spread 
around the mean value of log Al (hereafter we assume log 
to be base 10). Although real distributions are asymmet- 
ric, especially at higher luminosities, it is a useful first-order 
correction. Therefore, 



Pc,i(M) = G(logM;logMo,<,cT„,i), 



(5) 



where G{x\ fi, a) denotes a Gaussian distribution with mean 
/X and dispersion a. 

Following Mamdelbaum et al. (2004) the distribution 
for P^ci is assumed to reflect a double power law for the 
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number of satellites per halo: 



PncAM) 




M < Mo,i 

Mo,i < M < 3Mo,, 

otherwise 



Constants of proportionality arc chosen so that the 
is continuous and normalised and give the overall 
fraction ai = 0.2. Based on simulations we choose 
lowing canonical values for Mox- M(— 19) = 5x 10^^ 
M(-20) = 2 X 10^^/i"^Mo, M(-21) = 1 x 10"/i"^ 
fix av,i = 0.2. 

The total number of galaxies is given by the 
of the number of halos and the number of galaxies 
and therefore 

r?nc(M|Li) _ r/nc(M|Li) 



(6) 



function 
satellite 
the fol- 
h-^MQ, 
Mq and 

product 
per halo 



Pc{M\Li) = 



PNc(M|Li) 



/r?nc(M|Li)dM 

r]nNc{M\Li) 



Nc,i 

rjnc{M\Li) 

J ??nNc(M|Li)dM ~ N^ci 



(7) 

(8) 
(9) 



where Nc,i and A^NCi denote the total number density of 

galaxies belonging to a central/non-central class and belong- 
ing to the luminosity bin Li. The fraction of non-central 
galaxies Oi is therefore given by: 

Next we turn to the distribution of velocities. Let us 
then define the pairwise velocity distribution P{v,r\Li, Lj) 
as the probability that of all pairs between galaxies from 
luminosity bin Lj and galaxies from the luminosity bin Lj, 
we find one at separation r, whose velocity difference vectors 
has magnitude v. We emphasise that in our definition, this 
is an object-weighted quantity rather than volume weighted 
quantity. The pairwise velocity dispersion as a function of 
radius is then given by: 



2/ ir r ^ io P{v,r\Li, Lj)v dv 

a2{r\Li,Lj) = " . 

Jo P{v,r\Li,Lj)dv 

Similarly the correlation function is given by: 



(11) 



J-P{v,r\Li,LMv .^2. 

The pairwise velocity has two contributions. Single halo 
contribution comes from pairs residing in a single halo and 

dominates at small distances (< lh~^ Mpc), while the two 
halo contribution dominates at larger distances. 
Therefore we can write 

P{v, r, \Li,Lj) = P'\v, r\Li, \L, j) + P^''{v, r\Li, Lj). (13) 

We now investigate each contribution in turn. 



2.1 One-halo contribution 

The single halo contribution to the pairwise velocity distri- 
bution has contributions from pairs formed one from the 
central and one from the satellite galaxies (host-sat pairs, 
denoted as hs) and those coming from two satellite galaxies 
(sat-sat pairs, denoted as ss): 

P^''{v,r\Li,Lj) = P^,{v,r\Li,Lj) + P,,(v,r\Li,Lj). (14) 
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Figure 1. This figure sliows the predictions for tlic one-point 
velocity dispersion relative to the centre of the halo for satellites 
selected from different luminosity bins for the halo model. Three 
luminosity bins are plotted: —21 (solid line), —20 (dashed line), 
— 19 (dot dashed line). 



The host-sat contribution is proportional to 

J\s(w,r-|Lj,L,) oc y" [nc(M|Li)nNc(M|L,)-|- 

nc(M|L^)nNc(M|Lj)]Pi,(v|r,M)$i(r|M)77dM. (15) 

Note that even though formally we have two terms in- 
side square brackets of equation 15, if Lj < Lj first term is 
negligible by equation 6. 

The 1-particle velocity distribution P\^(v\r,M^ is as- 
sumed to be Maxwellian. This is true for isothermal spheres 
and is in reasonable agreement with simulations. Therefore 



P\v{v\r,Li) = M{v;K{Li)cTdin), 



(16) 



where M{x; a) denotes a Maxwell distribution correspond- 
ing to one-dimensional velocity variance aam, which we take 
to be 



: 90kms" 



M 



1012 ^-iMg 



1/3 



(17) 



For the sake of simplicity we neglect the velocity bias at 
this stage, simply setting by = 1, but we return to this point 
below when we discuss the simulations. Note that we can 
simply switch between calculating the distribution of one- 
dimensional components of velocity and the distribution of 
speeds (lengths of the 3-dimcnsional velocity vectors) by 
replacing the Maxwell distribution with the corresponding 
one-dimensional Gaussian. We will use the superscript ID 
to denote the distribution of one-dimensional components of 
velocity. 

Similarly, one can write the sat-sat contribution as: 



Pssiv,r\Li, Lj) oc J nNc(Af|Li)nNc(M|Lj) 

X P2v(y\M)<S>2{r\M)r]dM. 



(18) 



In this equation $2 denotes the probability that a pair of 
satellites is found at separation r. This corresponds to the 
self-convolution of the profile, which is most easily achieved 
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in the Fourier space: 

^2{r\M)(xr^J J dr^i{v\M)smc{kr) sinc(fcr)fc^dfc. 

(19) 

It should be emphasised that the halo occupations are 
stochastic and a probability distribution is required to fully 
describe their statistical properties. We have implicitly per- 
formed the appropriate averages in the above expressions. 
In particular, the number of pairs for a Poisson distributed 
variable p is given by {p{p — 1)) = (p)'^- 

The function P2v {v\AI) is used to denote the probability 
that the difference between velocity vectors hcis magnitude 
V. Because we assumed Maxwell distribution of velocities 
this is given simply by 



P2v{v\r,Li) = M{v;V2ad„ 



(20) 



The above equations can be combined in the following 
expression for the pairwise velocity distribution: 

[(i-„.)<.,p„(«it.)«,„,«i.,,P..,.|i,)+ 

(1 - aj)aiPG{M\Lj)P^G{M\Li)Pi^(v\Li)]^i{r\M)+ 
aiajPj,c{M\Li)P!,c{M\Lj)^2{r\M)P2v{v\Li,Lj)' 



dM 

V 



(21) 



This can be meaningfully compared to the one-point velocity 
distribution of galaxies at a distance r from the centre of its 
halo, that belong to the luminosity bin Li. In our notation, 
this is given by: 

Pi{v,r\Li) = J aiPNc(M|LO$i(r|M)Pi„(w|LOMdM. 

(22) 

Note that the contribution due to the peculiar motion of the 
entire halo is not present in the above distribution. 

Equations (21) and (22) show that true pairwise veloc- 
ities tend to favour higher velocities compared to the one- 
point velocity distribution for two reasons. First, in a halo of 
a given mass, the velocity difference vectors of sat-sat pairs 
have a higher dispersion (by a factor of V2 if we neglect 
velocity bias). Second, the number of pairs contributed by a 
single halo scales as the square of the number of its satellites 
and since heavier halos contain more satellites, the pairwise 
velocity dispersion puts more weight towards heavier halos. 
This is encoded in the weighting factor for true pairwise ve- 
locity distribution, which is r)~^dM compared to dM for the 
case of one-point velocity distribution. 



2.2 Two-halo contribution 

Two halo contribution comes from galaxies that are posi- 
tioned in different halos. For host-host pairs, the correla- 
tion function follows the correlation function of parent halos, 
while for host-sat pairs the correlation function would have 
to be convolved with the halo profile and for the sat-sat it 
would have to be convolved with the halo profile twice. How- 
ever, these effects are only important on small scales where 
one-halo term dominates. Hence, we ignore this effect and 
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Figure 2. This figure shows the results for the one-point velocity 
dispersion relative to the centre of the halo in the simulation. 
Lines are labelled as in Figure 1. 



assume that to a good approximation the three correlation 
functions are equivalent. 

The correlation function for halos follows that of the 
linear theory, with the complication that the more massive 
halos are more correlated. This is modelled in a simple man- 
ner by introducing a bias factor for halos of mass M, given 
by Sheth & Tormcn (1999) 



h{y) = 1 + 



2p 



(23) 



5c{l + v'p)' 

where u = [5c/cr{M)]'^ and u' — av. The spherical over- 
density at which a clump collapses Sc is 5c = 1.68 for 
Einstein-de Sitter model and fitted values of other parame- 
ters axe a = 0.73 and p = 0.15 (Mandelbaum et al. 2004). 
The cross-correlation between haloes of different masses is 
then given by: 



^i2(r|Mi,M2) = &„(r)6(Mi)6(M2). 



(24) 



The probability of finding a pair of halos of masses be- 
tween Afi,2 and Mi, 2 -|- dMi,2, separated by a distance be- 
tween r and r -|- dr is therefore given by 

2 



(1 + ^i2{r\Ml, M2)) 77(Mi)7?(M2)47rr 



(25) 



Differences between host-host, host-sat and sat-sat pairs 
arise when considering their velocities. The velocity differ- 
ence between host-host galax;y pairs arises purely due to 
velocity difference of their parent halos. The velocity dis- 
tribution for halos is often assumed Maxwellian (Sheth & 
Diaferio 2001). In our simulations, it was observed that to 
a good approximation the velocity dispersion is given by 
""halo = 455km s^'' , independent of halo mass. 

However, there are two factors that have to be taken 
into account when calculating the two-halo velocity disper- 
sion term. First, neighbouring velocities are heavily corre- 
lated and consequently the probability distributions for ve- 
locities of two halos cannot be assumed independent. Sec- 
ond, if we select halos with a nearby neighbour and com- 
pare its velocity dispersion to population of halos without a 
nearby neighbour we find that the former is higher. This is 
simply due to the fact that halos that have close neighbours 
are more likely to live in a denser environment and thus 
have higher velocity dispersion. These effects are difficult 
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to model and we combine them into an effective dispersion 
suppression factor s{r), so that the halo-halo probability for 
velocity difference can be written as: 

piUv\r, Ml, Ma) = M(v; sM^iahaio) (26) 
We use the following fitting formula for s(r): 

sir) = {i^y (27) 
1^ 1 otherwise, 

where ro = 50h~^ Mpc and v = 0.15. The relative bias of 
haloes of different mass is ignored in this approximation. 

When considering host-sat and sat-sat pairs we have to 
add the dispersion due to the satellite's peculiar motion with 
respect to its parent halo centre: 

pl'^{v\r, Ml, Ma) = M (v; [2aLioS{rf + cTdn,(M2)] 

(28) 

and 

pi''(i;|r,Mi,M2) = M{v; [2aL,,s(rf 

+ adm(Mi) + adn.(M2)]'/') . (29) 
Finally, we can assemble the 2-halo probability function: 

P'"'{v,r\Li,Lj) = J dMi J dMa 

X (^{1 - Qi)(l - aj)Pc{Mi\Li)PciM2\Lj)pli{v\Mi,M2) 

+ PNc{Mi\Li)Pc{M2\Lj)pt{v\M2, MuLj,Li) 
+ Pc (Ml ILOPnc (Ma \Lj )pt iv\Mu M^, Li, L^) 
+ QiQifkc(Mi|Li)fkc(M2|L,)pL^(v|Mi,M2)) 

X (1 + $i2(Mi, Ma))iViiyj47rr^ (30) 



3 BASIC PREDICTIONS AND COMPARISON 
WITH SIMULATIONS 

To check the validity of this model we compare it to a high- 
resolution numerical simulation. We use collision-less (dark 
matter only) simulation that is described in detail in Tasit- 
siomi et al. (2004), which has already been used for com- 
parison with the halo model (Mandelbaum et al. 2004). We 
used the largest available simulation with the box size of 
120h~^ Mpc and 512^ particles and the cosmological model 
parameters were the same as described in the previous para- 
graph. The resolution of the simulation limits the analysis 
to halos heavier than 10^^ Mq and thus we limit the dis- 
cussion to galaxies that are brighter than —19 in the r-band 
magnitude. 

The halos were identified using a variant of the Bound 
Density Maxima halo finder (Klypin ct al. 1999); the de- 
tails of this algorithm and the corresponding parameters can 
be found in Kravtsov et al. (2004). A catalog of dark mat- 
ter halos is constructed so that every dark matter halo can 
have zero or more subhalos. Every halo is assumed to host 
a central gala:xy. The maximum circular velocity is used as 
a proxy for halo mass for both halos and subhalos. Then 
the luminosities are assigned to each halo by matching the 



cumulative velocity function n{> Vmax) to the observed r- 
band cumulative luminosity function (Blanton et al. 2003). 
Since the mean SDSS redshift is ~ 0.1 the same redshift was 
used in simulation. More details can be found in Tasitsiomi 
et al. (2004). 

3.1 One-point velocity dispersion 

In Figure 1 we show the one-point velocity dispersion relar- 
tive to the centre of the halo. This is calculated using the 
equation (11), but using the one-point velocity distribution 
given by the equation (21). In other words it is the root 
mean square velocity of satellites at distance r from the 
centre of their halo. Note that the so-defined one-point ve- 
locity dispersion is different from host-sat contribution to 
the pairwise velocity distribution: the latter additionally re- 
quires the luminosity of the central galaxy to belong to a 
chosen bin. 

The velocity dispersion increases with increasing radius, 
even if we assume that the velocity dispersions is constant 
throughout a halo. This is simply a consequence of halo den- 
sity weighting. The one-point velocity dispersion shows the 
expected behaviour of increasing velocity with the increas- 
ing luminosity at small separations. For separations from 
the centre of the halo that are larger than cut-off radius of 
halos of mass Mo,i (for the luminosity bin Li), the velocity 
dispersion increases faster with radius. At large enough sep- 
arations all curves converge to the same one, following the 
distribution of velocities in halos weighted Mq, where only 
halos with large enough cut-off radius can contribute. 

Figure 2 shows the one-point velocity dispersion for 
galaxies in our simulation. We note that it does not show 
the expected behaviour. The more luminous galaxies live in 
heavier halos and thus they should be moving faster. This 
is a fairly robust prediction of the theory. The simulation 
shows that galaxies at the same distance from the centre 
move at an approximately constant speed at all radii, re- 
gardless of their luminosity. To understand this behaviour 
better we plot the one-point velocity dispersion relative to 
the centre of the halo in bins of the parent halo mass in Fig- 
ure 3. We split the entire galaxy catalogue into four logarith- 
mic bins according to the parent halo mass. We then split 
galaxies in each bin into three luminosity bins and plot the 
one-point velocity dispersion, where each individual point 
was plotted only if more than five galaxies contributed to it. 
We note the following: first, the approximation that the ve- 
locity dispersion is constant throughout the halo is roughly 
valid, in fact the velocity drops somewhat with distance as 
expected for an isotropic tracer in an NFW profile. This can 
also be a result of increased subhalo bias towards the centre 
(see e.g., Diemand et al. 2004). Second, the less luminous 
galaxies move faster than the more luminous galaxies even 
in halos of the same mass. This indicates a velocity bias that 
is increasing with decreasing luminosity. It is this bias that 
in Figure 2 conspires to nearly cancel the fact that one aver- 
age more luminous galaxies come from more massive halos. 

3.2 Pairwise dispersion 

We first investigate how individual components contribute 
to the pairwise velocity dispersion. We plot these in Figure 
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Figure 3. This figure shows the variation of the one-point velocity dispersion relative to the centre of the halo for galaxies belonging to a 
certain parent halo mass range. Line-styles are same as in the previous figures. The horizontal dotted line is the virial velocity dispersion 
of the mass corresponding the centre of the (log) mass bin. See text for discussion. Note that vertical axes have diflferent ranges. 



4. Thick lines are velocity dispersions of different compo- 
nents. The total velocity dispersion is the weighted combi- 
nation of these individual dispersions. To ease visualisation 
wo also plot individual components multiplied by the frac- 
tion of pairs of a given type. 

At separations larger than ~ 2/i~^ Mpc, the two-halo 
contribution dominates. At smaller scales there is no regime 
where one type would dominate. The (—21,-19) bin has 
a somewhat bigger contribution from the host-sat galaxy 
pairs, but not dominant enough to allow a clear separa- 
tion by distance. The sat-sat contribution has considerably 
higher dispersion than the host-sat contribution. These pairs 
come from massive halos that have many light satellites and 
since dispersion is a pair weighted quantity, these dominate. 

In Figure 5 we plot the combined contribution to the 
pairwise velocities as predicted by our model. At distances 
below lh~^ Mpc wc note that the galaxies selected from less 
luminous bins appear to move faster. This is explained by 
the fact that the fraction of host-sat pairs is decreasing as 
we go towards less luminous galaxies and sat-sat signal is 
almost independent of luminosity. This happens when the 
lower limit halo mass of satellites is low enough (i.e., the lu- 
minosity of galaxies involved is small enough) that the domi- 
nating sat-sat contribution, which probes M^^2{r\M)r) (see 
Equation (21)), approaches a constant. In Figure 6 we plot 



M^T] showing that sat-sat pairs effectively pick up masses 
around 10'^*/i~'^Mo (though this of course depends on the 
power-law index in equation (6)). The $2(f |M) factor addi- 
tionally biases towards higher masses with increasing r. 

When the two contributions are combined they result in 
a pronounced peak at around the typical size of the largest 
and therefore the most massive halos present. Such halos 
contain a lot of sat-sat pairs that move very fast and those 
dominate that region (at around 1 h ^ Mpc). At larger dis- 
tance we are dominated by many lighter halos with only a 
few members. 

In Figure 7 wc show the same total pairwise velocity 
dispersion as measured in our simulation. We note that the 
peak predicted by our simple model is also observed in the 
simulation. 

To summarise the conclusions from this section, we have 

found that the pairwise velocity dispersion consists of 3 com- 
ponents, host-sat, sat-sat and 2-halo pairs and that there is 
no regime where one component would clearly dominate over 
the others. 

The sat-sat component has the largest amplitude and 
drives the pairwise velocity dispersion peaks at around 
lh~^ Mpc. The amplitude at the peak is fairly constant 
and determined by the second moment of the mass func- 
tion, which is dominated by large clusters. Pairwise velocity 
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6v(-19) = 1.0, fev(-20) = 0.9, 6v(-21) = 0.6. Those values 
were chosen so that the model reproduced the correct be- 
haviour of the pairwise velocity dispersion. Results are plot- 
ted in Figure 8 (note that individual lines were shifted ver- 
tically by the same amount to make plot readable). The fit 
is remarkably good considering the simplicity of our model. 



3.4 EVaction of satellites 

A similar effect to the effect of velocity bias can be achieved 
by assuming a luminosity dependent fraction of satellites, 
i.e., a = a{L). As Figure 9 shows, increasing the satel- 
lite fraction leads to an increase in velocity dispersion, as 
expected. It also makes the peak at r = lh~^ Mpc more 
pronounced, since that peak is caused by satellite-satellite 
pairs. 

The numbers for the fraction of satellites in our simu- 
lation are 0.15, 0.20 and 0.23 for luminosity bins —19, —20 
and —21 respectively. Without introducing velocity bias it is 
possible to match the velocity dispersion of single-luminosity 
bins as shown in Figure 10. In this Figure we used the values 
as low as 0.06 for the most luminous galEixies, 0.17 for —20 
galaxies and 0.20 for faintest galaxies considered. 

However, as Figure 10 shows, the values of a that pro- 
duce a decent fit for galaxy pairs of a single luminosity bin 
considerably decrease the goodness of fit for cross-bins. We 
also note that decreasing the fraction of satellites decreases 
the height of the peak at around r ~ lh~^ Mpc. 



O I • • • ' • • • — 

0.1 1 10 

r [Mpc/h] 

Figure 4. In this figure we show dlfl^erent contributions to the 
pairwise velocity dispersion. The very thick solid line is the to- 
tal velocity dispersion, composed of host-sat pairs (thick solid), 
sat-sat pairs (thick dashed) and two-halo contribution (thick dot- 
dashed) . These three contributions are weighted by the respective 
number of pairs. The thin lines show velocity dispersion multi- 
plied by a fraction of a given component (so that thin lines add- 
up to the thick line). We plot the (-19,-19) (top) and (-21,-19) 
(bottom) pairs. See text for discussion. 

dispersion cannot be used to determine the typical halo mass 
of galaxies except on very small scales, where the number 
of pairs is small. This is true even for cross-correlation be- 
tween different luminosity samples, which in principle is a 
better way to select host-sat pairs: while this does increase 
the host-sat to sat-sat ratio we find that it still does not 
sufficiently suppress the number of sat-sat pairs to malce it 
a reliable tracer of host dark matter halo. Thus additional 
selections based on galaxy environment are needed if one 
wants to select a host-satellite sample only. Such additional 
selections may be difficult to model and reduce the size of 
the sample. 

3.3 Adding velocity bias 

We are now in position to rectify the main deficiency of 
our model with respect to the simulation, namely the veloc- 
ity dependent bias. This was achieved by using the follow- 
ing values of the velocity bias as defined in Equation (16): 



4 THE OBSERVABLE CORRELATION 
FUNCTION 

4.1 Correlation function and the halo model 

The pairwise velocity distribution we have been calculating 
so far is a basic distribution from which various other statis- 
tics can be calculated. It is not directly observable by itself. 
In the flat sky approximation, we see objects at various pro- 
jected radii and their redshift encodes both their distance 
and their velocity component. 

At large separations corresponding to two-halo contri- 
butions the galaxies experience infall that is well described 
by the linear theory beta model of Kaiser (1987). This re- 
sults in squashing of the correlation function along the line 
of sight. The amount of distortion is parametrised by a sin- 
gle parameter f3 = f /b, where / ~ il^^ is the linear growth 
factor and b is the linear bias. Ideally, one would like to in- 
corporate the infall by considering a distribution of infall 
velocities. Note that Scoccirnarro (2004) has shown that no 
distribution of infall velocities reproduces the beta model. 
Nevertheless, since beta model appears to be a good phe- 
nomenological model for description of infall, we include it 
here. 

First we rewrite the probability of finding an object 
between r and r -|- dr and having velocity v between v and 
V + d.v into cylindrical coordinates rproj and z: 

Pcyiir.r.i,z,v) = p{^^l~T7^v)^^^^^^ (31) 
The distribution in the rproj — tt plane is then given by 
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Figure 5. This figure sliows tiie total prediction for tile pairwise 
veloeity dispersion for various luminosity bins shown in the leg- 
end. Distances below \ Mpc are dominated by the one-halo 
term, while the two halo term dominates at larger distances. 



-Pobs( ' p 



■Pcyl (»'p 



,«,i?o(7r - z)) 



(32) 



We then calculate the observable correlation function 

by summing up contributions from 1-lialo and 2-halo terms, 
following a procedure similar to many other works (White 
2001; Seljak 2001; Kang et al. 2002; Cooray 2004): 



proj,7r) -|-^^''(rp 



(33) 



The one-halo contribution is calculated from the Equa- 
tion (32). This implicitly assumes that these systems are 
completely virialised and therefore not affected by the infall. 
The two-halo contribution, on the other hand, is calculated 
using standard dispersion model, which takes the following 
simple form in Fourier space: 



P^(k) = P^{k) 



(1 + fip? 



(34) 



Here and denote the power spectrum in redshift 
space and real space respectively and /i is angle to the line 
of sight. The /3 parameter is calculated using bias predicted 
for a given luminosity bin from our model and velocity dis- 
persion is given by the equation (11). 




Figure 6. This figure shows r;M^ as a function of M and shows 
which mass range is picked up by the sat-sat pairs in our model. 
See text for discussion. 
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Figure 7. This figure shows the total pairwise velocity dispersion 
as measured in the numerical simulation. 



4.2 Comparison to observations 

Jing & Borner (2004) have found that less luminous galax- 
ies in the 2dF redshift survey catalog appear to move faster. 
This conclusion is reached as follows. The correlation func- 
tion of galaxies in the 2dF catalog is calculated for galaxies 
selected from various luminosity bins. This correlation func- 
tion is used to calculate the redshift-space power spectrum 



Pairwise velocities in the Halo Model: Luminosity and Scale Dependence 9 




(-21,-21) 
(-20,-20) 
(-19,-19) 



r [Mpc/h] 




(-21.-20) 
(-20,-19) 
(-21,-19) 




(-21,-21) 
(-20,-20) 
(-19,-19) 



r [Mpc/'h] 




(-21,-20) 
(-20,-19) 
(-21,-19) 



0.1 



10 



0.1 



10 



r [Mpc/h] 

Figure 8. This figure shows the total pairwise velocity disper- 
sion as measured in the numerical simulation compared to our 
model when velocity dependent bias has been added to the model. 
Jagged thin line are simulation data and smooth thick line the 
model prediction. Lines have been shifted vertically to ease visu- 
alisation. 



0.1 



1 



[Mpc/h] 



Figure 9. In this figure we plot the eflfect of changing the satel- 
lite fraction on the velocity dispersion. All lines correspond to 
(-19,-19) pairs the satellite fraction a set to 0.4 (dot-dashed), 

0.2 (dashed) and 0.1 (solid). Increasing the a parameter increases 
the overall velocity dispersion, but also makes the lh~^ Mpc 
bump more pronounced. 



r [Mpc/h] 

Figure 10. Same as Figure 8, but with variable fraction of satel- 
lites used to match single-luminosity bins, instead of luminosity 
dependent velocity bias. See text for discussion. 



P'^(k). For each value of k the data is fit by fixing (3 = 0.45 
and fitting for ai2 and P{k). 

There are three features in our model that can lead to 
an explanation of this effect. First, it is possible that the 
effect of pair weighting gives enough biasing towards more 
massive systems to explain the effect. As shown in Figure 11 
the recovered velocity dispersion is non-monotonic, allowing 
for the possibility that velocity dispersion decreases with lu- 
minosity over some range, although the effect is not very 
pronounced and depends on the scale. Second, the addition 
of the luminosity dependent velocity bias works further in 
making less luminous galaxies appear to move faster. And 
third, allowing for satellite fraction to decrease with lumi- 
nosity also leads to a decrease of velocity dispersion. 

In order to meaningfully compare our predictions with 
Jing & Borner result we mimic their procedure. We predict a 
rcdshift space correlation function using our model, calculate 
the power spectrum and then fit the dispersion model to it. 
We also follow their procedure by fixing the (3 parameter to 
0.45. 

The main results of this exercise are plotted in Figures 
11 and 12. We plot the recovered value of the ar2 parameter 
as a function of the wave vector k for three different lumi- 
nosity bins. Figure 11 shows the behaviour for the simple 
model discussed throughout this paper. We note that the 
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Figure 11. This figure shows the recovered (T12 pajrameter follow- 
ing the prescription used to analyse the data in (Jing & Borner 
2004). Line styles as before. 



5 DISCUSSION AND CONCLUSIONS 

In this paper wc have investigated the pairwise velocity dis- 
persion in the halo model and compared it to simulations. 
The distribution of velocities for galaxies separated by a cer- 
tain distance in the real space is a pair weighted quantity and 
thus does not reflect the velocity structure of halos typically 
associated with galaxies of a given luminosity. In particular, 
for the pairs selected from either unequal luminosity bins 
or from equal luminosity bins (e.g., (—20,-19), (-19,-19), 
etc.) the main contribution comes from the sat-sat pairs. 
Since the halo mass distribution of satellites follows the halo 
mass function, which is exponentially suppressed on high 
mass end and has a power law at the low mass end, it probes 
a reasonably restricted range of parent halo mass around 
10^*h~^ Mq and is nearly a constant. This, in combination 
with the fact that if one has a larger number of galaxies the 
number of sat-sat pairs relative to host-sat pairs increases, 
results in velocity dispersion which is not very luminosity 
dependent. The pairs selected so that one galaxy is more 
luminous that the other are also dominated by sat-sat pairs. 




0.2 



0.5 



kproj [h/Mphc] 



Figure 12. Same as figure 11, but a luminosity dependent veloc- 
ity dispersion has been added to the model. 



plot looks qualitatively the same as Figure 5 with an in- 
verted horizontal axis. Here too, there is a typical bump at 
k ~ Ift/Mpc. The pair-weighting can result in less luminous 
galaxies appearing to move faster, but only at fc > 2/i/Mpc 
or so and tlicrforc it is not clear whether this effect alone 
provides satisfactory explanation. 

In Figure 12 we repeated the analysis, but added a lu- 
minosity dependent velocity bias that was used in section 
3.3. Velocity biasing has the expected effect on the inferred 
velocity dispersions: galaxies of lower luminosity can indeed 
have a lower ai2 even at wave vectors as large as 1/i/Mpc. 

We note that a similar analysis of the SDSS data (Li 
et al. 2005) shows qualitatively similar behaviour to our 
predictions plotted in Figure 12: a wide bump at around 
1 — 2ft/Mpc that becomes less pronounced for higher lu- 
minosities. The most luminous bins in that paper show an 
increase of ai2 with increasing k. We do not observe this 
behaviour here. 



We have made predictions for redshift space correlation 
function of galaxies and fit it with the dispersion model, 
mimicking the observation procedure employed in Jing & 
Borner (2004) and in Li et al. (2005). In the latter work 
the measured ai2 in SDSS data has a qualitatively similar 
dependence on the wavenumber as predicted by our model. 
One of the main conclusions of these papers is that velocity 
dispersion declines with luminosity, which appears to con- 
tradict simple halo models. As mentioned even the simplest 
models do not suggest that there is a strong luminosity de- 
pendence, since most of the contribution comes from sat-sat 
pairs. Two other effects can give rise to velocity dispersion 
decreasing with luminosity. One is luminosity dependent ve- 
locity bias. In the simulation the effect that more luminous 
galaxies appear to move slower is clearly evident: in the halos 
of the same mass, the more luminous subhalos tend to move 
slower than their less luminous counterparts. Tidal strip- 
ping and tidal disruption have been suggested to explain 
the anti-bias of subhalo distribution (see e.g., Kravtsov & 
Klypin 1999; Diemand et al. 2004; Zentner et al. 2005) to the 
dark-matter distribution in a given halo and to explain the 
velocity bias of the subhalo population. Subhalos with lower 
orbital energy do not survive and one is left with halos that 
are on average faster. We note however that in these dark 
matter simulations one uses a simple prescription to link 
the luminosity of a galaxy with the circular velocity of the 
subhalo, which may be more complicated in the real world. 
Second effect that can give rise to velocity dispersion de- 
creasing with luminosity is if satellite fraction also decreases 
with it. There is already some observational evidence for this 
from halo modelling of luminosity dependent galaxy corre- 
lation function Zchavi et al. (2005) and this explanation has 
also been put forth by Jing & Borner (2004). Both of these 
effects can therefore explain the observed trends and can be 
included in halo models. In summary, using auto and cross 
correlations between luminosity bins in redshift space to ex- 
tract PVD is unlikely to provide us with simple information 
such as the mass and mass profile of the central galaxy halos, 
but is instead telling us more about the fraction of galaxies 
that are satellites and their velocity bias inside clusters. 
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